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Abstract 

We study the perturbative unitarity bound given by dimension six derivative interactions consisting of 
Higgs doublets. These operators emerge from kinetic terms of composite Higgs models or integrating out 
heavy particles that interact with Higgs doublets. They lead to new phenomena beyond the Standard Model. 

One of characteristic contributions by derivative interactions appear in vector boson scattering processes. 
Longitudinal modes of massive vector bosons can be regarded as Nambu Goldstone bosons eaten by each 
vector field with the equivalence theorem. Since their effects become larger and larger as the collision energy 
of vector bosons increases, vector boson scattering processes become important in a high energy region around 
the TeV scale. On the other hand, in such a high energy region, we have to take the unitarity of amplitudes 
into account. 

We have obtained the unitarity condition in terms of the parameter included in the effective Lagrangian 
for one Higgs doublet models. Applying it to some of models, we have found that contributions of derivative 
interactions are not so large enough to clearly discriminate them from the Standard Model ones. We also study 
it in two Higgs doublet models. Because they are too complex to obtain the bound in the general effective 
Lagrangian, we have calculated it in explicit models. These analyses tell us highly model dependence of the 
perturbative unitarity bounds. 
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1 Introduction 

Almost all experimental results are consistent to what the Standard Model (SM) predicts. The recent obser- 
vation of a new boson [1] have been also predicted by the model as the Higgs boson. Measured properties of 
the observed particle are still consistent to predicted ones, so that the SM is established by the observation. 
However, the model seems to be the effective theory describing physics below the electroweak (EW) scale 
since it includes several theoretical and experimental problems maybe explained at the TeV scale. 

One of those mysteries is how to break the EW gauge symmetry. Since the symmetry is broken by hand 
in the SM, the model tells us what happens in experiments up to 0(100) GeV but does not tell us why and 
how the EW symmetry is broken. Therefore, we expect the SM should be extended to a model including other 
sectors responsible for physics of the electroweak symmetry breaking (EWSB) and other veiled phenomena. 
If the scale is much higher than the EW scale, the observed Higgs mass requires subtle mechanisms or fine 
tuning. Therefore we expect something new appears at several TeV region. Eventually, new structure is 
expected to appear at the TeV scale in order to obtain the EW scale without those artificial constructions. 

Many of models describing physics beyond the SM have been proposed. Higgs sectors of these models 
are extended, and they break the EW symmetry with a certain mechanism. Effects of these new sectors could 
be firstly observed as deviations from the SM ones via higher dimensional operators at the scale lower than 
their original scales. 

When we consider physics beyond the SM with extended Higgs sector, its low energy effective theory 
probably include dimension six derivative interactions as a part of higher dimensional operators. These op- 
erators have two origins: expansion of their kinetic terms if the Higgs doublet is realized as a part of pseudo 
Nambu Goldstone (NG) field; integrating out of heavy new scalar/vector bosons that interact with the Higgs 
field. The latter case appears even in models including elementary Higgs. 

If the Higgs boson were removed from the SM, the Higgs sector would be described by SU (2)l x 
SU (2)r/SU (2)y nonlinear sigma model. Derivative interactions of NG fields emerge from the kinetic term. 
These interactions contribute to longitudinal mode scatterings of massive gauge bosons with the equivalence 
theorem and cross sections of these processes become larger and larger as energy increases. They finally be- 
come so large as to violate the perturbative unitarity about 1 TeV [2]. Of course, the recent observation of the 
Higgs boson told us absence of the unitarity violation and reliability of the SM description even much above 
the TeV scale. We confront the similar problem in studying derivative interactions of Higgs doublets. In this 
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paper, we clarify scales where the given perturbative description is available with derivative interactions of 
the Higgs doublet in several models. 

The rest of this paper is organized as follows. In Sec. 2, we study the unitarity bound given by derivative 
interactions in one Higgs doublet models (lHDMs). It is extended to the case of the two Higgs doublet 
models (2HDMs) in Sec. 3. In both of these sections, the unitarity violation scales are explicitly calculated 
with several models. Finally, our study is concluded in Sec. 4. 



2 Unitarity of derivative interactions on one Higgs doublet models 

The perturbative unitarity bounds given by the derivative interaction are discussed on lHDMs. Firstly, we 
derive the formula of the unitarity bound and investigate its general properties. Then, results are applied to 
explicit models. The formulae of the perturbative unitarity are shown in App. A. 



2.1 Formulae and general properties of the unitarity bound 

The effective Lagrangian of derivative interactions in 1HDM is 3 

H T 

* ^Yf2 d ^ H ) d ^ H ) + Yf2 ^ (H^H) , (2.1) 

where / is a scale related to new physics. The second operator should include covariant derivatives. However, 
gauge fields introduced by this kind of operator are not studied in this paper so that covariant derivatives are 
replaced with partial derivatives. Since the latter term violates the custodial symmetry, our analysis is based 
on the Lagrangian with c T = 0. 

Since we consider only four point scatterings given by Eq. (2.1), the vacuum expectation value of the 
Higgs boson play no role in the following calculation. Therefore, we use 

H = ^ , + = (C- , (2.2) 

where C + /N is a charged/complex-neutral scalar field. The charged scalar and imaginary part of the neutral 
scalar are respectively eaten by W ± and Z bosons. Using the above notation, the following amplitudes are 
obtained 4 : 

J?{C + C- -»• C+C~) =^(AW t AW 1 ) 

=^V, (2-3) 

J?(C + C- -»• AW 1 ) =4tC H , (2.4) 
/ 

where s and t are the Mandelstam variables and particles are considered as massless, i.e., s + t + u = 0. 

Following Ref. [2], we construct matrices with partial wave amplitudes. The largest eigenvalue of these 
matrices give us the strongest bound to the perturbative unitarity. We have found zeroth modes produce the 
strongest bound in 1HDM, so that we focus on the case. With the formulae in App. A, the bound is given by 
the largest eigenvalue of the following matrix: 

(M Q {C + C~ ^C + C-) M (C + C~ -+NN*)\ _J_(c h /2 c H \ 

V M (NNl -> C+C-) M {NN^ -»• AW+) J ~ I6n f 2 \ c H c H /2 J ' ( } 



namely, the bound is 



J^^- (2-6) 



3 Using the field redefinition H —*H+ {a/ f 2 ){H^H)H, where a is chosen as an appropriate value, any other dimension six 
derivative interaction of the Higgs doublet can be expressed with these kinds of operators given here [3]. 

4 The EWSB produces corrections of 0(v 2 // 2 ). However we neglect them as neglecting 0(m 2 y Z /E 2 ) corrections associated with 
the equivalence theorem and dimension eight operators because they are subleading effects in this study. 
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2 
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2/3 
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W+W+ - 




1 


1/2 



Table 1: Cross sections for vector boson scattering (VBS) processes in the unit of a(W^W L — > hh). In the 
column of Full/Central, cross sections of VBS sub processes with/without the central region cut are shown. 



Assuming that derivative interactions are purely given by kinetic term of the nonlinear sigma model, the 
conservative cut off scale is expressed in terms of the decay constant, i.e. A ~ 4nf 5 . Using the relation, the 
unitarity bound is related to the cut off scale as 

' 1 (2-7) 



A 2 3kc h ' 
Therefore, if the relation 

c»<± (2.8, 

is satisfied, models reach the cut off scale before accessing the unitarity violation scale. Then, the effective 
Lagrangian, Eq. (2.1), is available up to the cut off scale. On the other hand, if the coefficient, c H , is much 
larger than unity, the unitarity violation scale becomes comparable to the scale of new physics, /, so that the 
description of the effective Lagrangian is invalid even in the energy region around /. In the case c H is 0(1), 
the unitarity violation scale exists between the scale of new physics and the cut off. Most of examples shown 
later are involved in this case. Around the unitarity bound, we have to include resonance effects, see, e.g., 
Ref [5]. It is therefore necessary to clarify valid energy scales for the description on each model. 

We apply the result to cross sections of the Higgs boson and longitudinal modes of massive gauge bosons 
with the equivalence theorem. Since these scatterings are dominated by the coefficient, c H , with the custodial 
symmetry, all of cross sections are proportional to each other. Here, we focus on only the process W^W£ — > 
hh, and relations with the others are shown in Tab. 1 . Considering this sort of processes, the importance of 
the central region 6 has been pointed out in Ref. [6], so that we also show ratios between the cross section of 
the Higgs pair production and those of the other processes with the central region cut. The cross section of 

w l w l -> hh is 

a(W+Wr -> hh) =— ( C) Z — * L ? Xl ° 6 [fb]. (2.9) 
v L L ' 3271 \ f V 9s s[TeV} 2 1 J 

For the process, Fig. 2.1 shows the region where the perturbative unitarity is violated . 

Assuming that cross sections reach the above bound at y/S = 3 TeV, the decay constant is 



' ; 733 [GeV]. (2.10) 



167T 

If c H ~ 1 , the effect of the derivative interaction in the process becomes comparable with the SM one about 
y/S = 2 TeV, where the cross section is 3 x 10 4 fb without the central region cut, see Ref. [6]. The value of / 
is typically related to new particle masses. For example, in the little Higgs scenario [7], the top partner mass 
is given by 0(f). From the viewpoint of the fine tuning, / is required to be below about 1 TeV. 



5 If UV completions are specified, the generalized dimensional analysis may produces lower cut off scales [4]. 
6 This region is defined as cos 6 e [— 1/2, 1/2] in detectors, where is an angle from the beam axis. 



3 




500 1000 1500 2000 2500 3000 
y/i [GeV] 



Figure 2.1: The upper bound of the cross section for W^W^ hh with the perturbative unitarity condition. 
The horizontal axis is the collision energy of this VBS sub process. In the upper shaded region, the unitarity 
is broken down. The black, dark gray and light gray lines are the cross sections where f/Vc^ = 500,750 and 
1000 GeV, respectively. For the other processes, the bounds can be obtained with the shift of the vertical axis 
by factors given in Tab. 1 . 



2.2 Examples with explicit models 

In the rest of this section, we study the unitarity bounds on two models: the minimal composite Higgs 
model [8]; the littlest Higgs model with T-parity [9]. The latter model is studied in Refs. [10, 11]. In 
Ref. [11], several Little Higgs models are also investigated 7 . Since the normalization of decay constants 
can be changed, the combination f 2 /c H is meaningful. We here follow the normalization given in original 
papers. Decay constants have physical meanings through masses of additional massive vector bosons and 
fermions in each model. 

The Higgs doublet is embedded so as to reflect the custodial symmetry. Therefore the operator violating 
the symmetry does not appear in both models. 



2.2.1 The minimal composite Higgs model 

This model is described by SO(5)/SO(4) nonlinear sigma model including four NG fields. They are identified 
as the Higgs doublet. 
The Lagrangian is 

if = ^(dZ) Jf (dZ), (2.11) 

with 

Mh/f]h/h 

cos[h/f] J 

where h is the real scalar multiplet of four NG bosons and h is its norm. Expanding these trigonometric 
functions, it is obtained that 

c H = \. (2.13) 



7 They obtained the unitarity bounds with all NG bosons. On the other hand, we focus on the Higgs doublets because other NG 
bosons are too heavy to treat as massless particles. Our results are conservative compared to theirs. 
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With the Eq. (2.6), the relation between the decay constant and the scale of the unitarity violation is 



s 



\6n 



(2.14) 



f 2 
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Assuming that the perturbative unitarity is violated at 3 TeV, the decay constant is about 750 GeV . On the 
other hand, if the decay constant is chosen as 500 GeV, the perturbativity is preserved up to about 2 TeV, 
where the cross section of W^W[ — > hh is 7 x 10 5 fb. In the case, the cross section of the Higgs boson pair 
production is one order of magnitude larger than that given by the SM. However, it is challenging to observe 
this process because the main decay mode is into four bottom state suffered from QCD background. 

2.2.2 The littlest Higgs model with T-parity 

Derivative interactions on the littlest Higgs model with T-parity is shown below. The scalar fields are described 
by SU (5) /SO (5) nonlinear sigma model which include 14 NG bosons. 
The kinetic term of this model is 



We omitted the other NG bosons since they don't contribute to the current analysis. Extracting the derivative 
interaction from the kinetic term, we obtain 




This result is consistent with previous works [10, 11]. 

If we suppose that / = 750 GeV, the perturbative unitarity is preserved up to about 4 TeV. Hence this 
model description is valid in higher energy scale while signals of the derivative interaction are smaller than 
the previous model. For / = 750 GeV, the cross section of W^W[ — > hh in this model is almost corresponding 
to the line where f/Vc 11 is 1000 GeV. 

3 Unitarity of derivative interactions on two Higgs doublet models 

In this section we extend the previous discussion to dimension six derivative interactions including two Higgs 
doublets. 

The modification is straightforward, and the prescription is also simple. However, formulae become too 
complex because of many degrees of freedom (DOF). Then we cannot obtain the formula of the strongest 
bound like Eq. (2.6) with the largest eigenvalue of a matrix that consists of partial wave amplitudes. Since the 
matrix can be diagonalized in individual models, three models are shown as examples. 

In this section, we consider processes whose initial states are electromagnetic neutral. Matrices producing 
the unitarity bounds for singly or doubly charged initial states are also shown in App. C. 




(2.15) 



where J7 is NG field. The Higgs doublet is assigned in the NG field as 




(2.16) 
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3.1 Formulae and general properties of the unitarity bound 

The analyses in this section are based on the following effective Lagrangian: 

H H 
q>~. C UU n H , c U\2 (n H , n H \ 
J_ ^2 _C 'llll + _ ^2 _ V C 'lll2 + C 'll2lJ 

i c U22 r> H i C 1221 rfl i C 1212 ( n H , n H \ 



, C 2221 , ,>// , ,>// \ , C 222 _' n H 
TT \ U 222\ C/ 2212J T? - ^2222 



22: 

y2 V^2221 1 ^2212; 1 y2 



where 



+ ^im + ^i0[ 221 + ^(0[ 212 + T 2m ), (3.1) 



°V« = \ + 8 ik 8 jl d{H i Hj)d{H * Hl) ' (3 - 2) 
°V« = i-u^ ^HDiflffm), (3.3) 



and 



= (^) ' ^ = ^ ^) • (3 " 4) 
To impose the custodial symmetry, the above coefficients are real and follow relations derived in App. B: 

3cfi22 + c 1221 ~~ c 1212 =0' (3-5) 
c fl22 + c 1221 + c r212 =0- (3-6) 

The 2HDMs require mixing angles to get mass eigenstates of scalar fields. In this paper, we use the 
equivalence theorem and focus on only derivative interactions, that is, masses of scalar fields are neglected. 
In this case, the perturbative unitarity bound is independent of mixing angles. This is also true for models 
including N Higgs doublets. 

Finally, the unitarity bound is 

A<J!L_ (3 7) 

J I*-' max I 

where C max is the largest eigenvalue of the matrices given in App. C 8 . 

As we will see later, the largest eigenvalue |C max | can be as large as about ten. In this case, the unitarity 
bound becomes quite strong and leads to an interesting remark. Suppose that, for instance, the pair production 
of a heavy particle whose mass is 0(f) in VBS processes, the energy scale where the pair is produced could 
be as large as the unitarity violation scale. This means we couldn't discuss the kind of process with these low 
energy descriptions. 



3.2 Examples with explicit models 

We study consequences of the above result with several models including two Higgs doublets. Following 
three models are studied: the bestest little Higgs model [12]; the UV friendly little Higgs model [13]; an 
inert doublet model. The first and second ones are composite Higgs models and the last one is a toy model 
including elementary Higgs doublets. 

8 Eq. (3.7) and the viewpoint explained below are also stated in Ref. [11]. 
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3.2.1 The bestest little Higgs model 



The bestest little Higgs model is a little Higgs model which includes two Higgs doublets. Scalar fields are 
implemented as SO(6) x SO(6) / SO(6) nonlinear sigma model including 15 NG bosons. The normalization 
of the kinetic term is the same as Eq. (2.15), and the NG field is 

hi /i 2 N 

n = —\-h T 1 I , (3.8) 



V2 



-hi 



where h\^ are real scalar multiplets considered two Higgs doublets and the other NG bosons are eliminated. 
In this model, Higgs doublets interact with heavy gauge bosons and a singlet scalar. The masses of heavy 
gauge bosons depend on the other decay constant that is larger than / in order to avoid constraints from the 
electroweak precision measurement (EWPM). Thus effects coming from the heavy gauge bosons are tiny, 
then we neglect them. The interaction with a singlet is required to obtain a collective quartic coupling. For 
simplicity, we introduce the following terms to see the effect: 

2 

=Sf ff D -^<7 2 + A/<7(# 1 t # 2 + H.c), (3.9) 
where a is a neutral singlet scalar 9 . Including this contribution, coefficients of derivative interactions are 

<fm=^ <f 112 =0, (3.10) 



2' 



1 « 1 



cf 122 — 0, c ?22i ~4 + c(7 ' c i2i2 — 4+ cC7 > (3-11) 

c 2221 =0, <?2222 (3-12) 



2' 
1 

4' 4' 



C U22 — 0, cf 2 2i — T 5 c r212 — ~~ 7' (3.13) 



where 



1 2 r 4 

c a = ±4-. (3.14) 

The unitarity bound depends on the value of c a because the largest eigenvalue is a function of it. For < c a < 
1/8, the bound is 

(315) 

At c a = 0, it is bounded as 

s 



, 2 <8tt, (3.16) 



and it becomes weak as c a increases. At c a = 1 /8, the bound is the weakest: 

J_ < 16; 
f 2 - 15 



i<^- (3-17) 



In the region, 1/8 < c a , the bound is 



where the right hand side decrease as c a increases and the bound becomes the same as the case of c a = at 
c a = l/l. 



9 In the original paper [12], m a = \/Ag5 + A56 / and X = ^ 65 ^ 56 ■ 
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The unitarity bounds of W^W L — > hh and — > W^W^ are displayed below. We define mass eigen- 
states, h and W^, as follows: 

h 

— =Af cosa+Af sina, (3.19) 
V2 

W£ =Cf cos J8 + cf sin J8 , (3.20) 
where Nf is the real part of N, and a and j3 are mixing angles. Unitarity bounds for these processes are 

2 



o(W+W L -+ hh) =J— A B h (^f £ ^^4^, (3.21) 



max 



a(W+W+ =^B w (Pf < ^^f, (3-22) 



-max 



where 

fi*(a,/5) =J (1 + (1 +2c ff (l -c 4p ))c 2{a _ p) +2c a s 4p s 2{a _^) , (3.23) 
B w (l5)=^(\+c a (\-c 4 p)). (3.24) 

Parameters c x and s x are cos^: and sinx, and C max = (2 — c a )/2 for < c CT < 1 /8 and C max = (1 + 7c CT ) /2 for 
1/8 < c° '. If a = j3 is satisfied, it is called the decoupling limit, and we get the relation: Bh{fi,fi) = S w (j3). 

The perturbative unitarity bounds of W^W[ — > hh are shown in Fig. 3.1. In order to see the effects of 
new parameters, we fix the decay constant as 750 GeV. The shaded regions in these figures are changed by the 
mixing angles because the cross section depends on the angles. However, the unitarity bound itself depends on 
only the coefficient, c° ' . Hence we can see that the energy scales where the cross sections intersect the unitarity 
violation regions are independent on the angles, e.g., y/S ~ 1.9 TeV for c a = 1. For j3 = and a — j3 = it/ 6, 
the cross sections don't depend on the value of c a \ thus, we have only one line but still intersecting points are 
the same. 



3.2.2 The UV friendly little Higgs model 

The UV friendly little Higgs model also includes two Higgs doublets in a part of 14 NG bosons given by 
SU(6)/Sp(6) nonlinear sigma model. The normalization of the kinetic term is also the same as Eq. (2.15). 
Since we study only Higgs doublets, NG field 17 can be considered as follows 10 : 



n = - 



e{H[-Hl 



V 



-e{H 1 -H 2 ) H 1 +H 2 



-HI -HI 



-b(h;-hz) 



e{Hl-Hl) 



\ 



J 



(3.25) 



where e is the totally antisymmetric tensor, £ 12 = 1. Since contributions given by heavy new particles can 
be ignored in this model, only derivative interactions generated by the kinetic term are taken into account. 
Following coefficients of derivative interactions are produced by the kinetic term: 



c im 

c 1122 
c 2221 
c 1122 



=4, 

=1, 
=0, 

=1, 



c 1112 
c 1221 
c 2222 
c 1221 



=0, 
=4, 
=0, 



c 1212 ' 



'1212 



= 1. 



(3.26) 
(3.27) 
(3.28) 
(3.29) 



Assignment and normalization of NG bosons given by the original paper are different from ordinary prescription of the nonlinear 
sigma model. Since we study only the part of two Higgs doublets, normalization of these field are changed in order to study with the 
canonical normalization. 
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Figure 3.1: The perturbative unitarity bounds of W^W£ — > hh for various c a and mixing angles. The decay 
constant is fixed as 750 GeV. The Mixing angles are set as (j8,a — j8) = {it/ '6,0) (upper left), (0,7i/6) (upper 
right), (7i/4,7i/6) (lower left) and (tt/6, 71 / 4) (lower right). The light gray, dark gray and black lines are cross 
sections for c a = 1,0 and 1/8, respectively. Unitarity violation regions depend on the value of c c , and their 
brightness correspond to each line. 
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Figure 3.2: The perturbative unitarity bounds of W^W£ — > hh in the UV friendly little Higgs model. The 
horizontal axis is the collision energy of this VBS sub process. In the upper shaded region, the unitarity is 
broken down. The black, dark gray and light gray lines are the cross sections corresponding to / = 500,750 
and 1000 GeV, respectively. 



These coefficients apparently violate the custodial invariant conditions, Eqs. (3.5) and (3.6). However, tree 
level contributions to p parameter cannot occur because of Z2 symmetry. With these coefficients, the strongest 
bound is 

jz<n. (3.30) 

Assuming that the perturbative unitarity is violated at 3 TeV, the decay constant, /, is determined as 1 .7 TeV. 
This value looks large in the viewpoint of the fine tuning as we have already seen. On the other hand, if the 
decay constant is about 1 TeV, the unitarity is broken below about 1 .7 TeV. 
The unitarity bounds of W£W£ -)■ hh and W L + W L + -> W L + are 

o(W+W L ^hh)=JL;<^ (3.31) 

<j(W+W+ -»> W£W?) =^jj < j f (3.32) 

Note that the cross sections have no mixing angle dependences because of Z2 structure; only one Higgs doublet 
gets the vacuum expectation value. These bounds to cross sections are shown in Fig. 3.2. It corresponds to 
the case c° = 15/7 for the bestest little Higgs model. The unitarity bound of this model is severe because the 
largest eigenvalue is quite larger than the previous models. 
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3.2.3 Inert doublet models with odd scalars 

We investigate the following Lagrangian consisting of elementary scalar and vector fields: 

2 

se d - $ + Ao/^o (h\h 2 + H.c. 

„2 



- ^ mi + Wl (H\o a H 2 + H.c.) 

- m 2 sL (j>f(j>^ + V2X L f ((j>f(Hfo 2 (J a H 2 ) + H.c.) 

+ ^V • V + g Vb • (mffH 2 + H.c.) 

+ m^yj • V s + v/2g (iV s f • Hfo 2 fH 2 + H.c.) 

2 

+ ^VZ-Vl + gLVt- (iH\c a d # 2 + H.c.) . (3.33) 



Scalar fields, 0o, 0£ and are respectively lo, 3o and 3i representations of SU(2)l x U(1)y, and vector 
fields, Vb, Vs and V£ ! are respectively lo, li and 3o representations of the gauge symmetry. We suppose that 
these new particles and H 2 are odd under an additional Z 2 symmetry, and Hi and the other SM particles are 
even under the discrete symmetry These choices of couplings and masses for and <j>f, and Vb and Vs are 
required to respect SO(4) symmetry 11 . This set up suppresses contributions to oblique corrections. 
Integrating out heavy particles, following coefficients of derivative interactions are obtained: 



cfm =0, cf 112 =0, (3.34) 

C H22 = - 2s l + 3v + 3v L , cf 221 =s - 2s L + 3v , cf 212 =sq - 2s L + 3v L , (3.35) 

c 222l =0, cf 222 =0, (3.36) 

c fi22 = - v o + v L , c\ 22X =-s l - v l , c[ 2l2 =s L + v , (3.37) 



where 



\m v0 J \m vL J 

Even if additional particle exist, their contributions are included into these four coefficients. We cannot dis- 
criminate these multiple contributions from a large contribution given by a particle. Using these coefficients, 
eigenvalues of the matrix (C.l) are fortunately obtained as the following simple forms: 





so - 7sl + 3v + 3v L 
2 


(3.40) 


4 


7s -9s L + 9v + 9v L 
2 


(3.41) 


4 


so ~ 3s L + 3v - 9v L 


(3.42) 


2 




s Q - 3s L - 9vo + 3v L 


(3.43) 


2 


4 ± 


s + s L + 3v + 3v L 
2 


(3.44) 


cL 


sq + 9s L - 9v Q - 9v L 
2 


(3.45) 



"This structure should be broken by renormalization group running of them even if the UV completion possesses the structure. 
We assume that this 50(4) symmetry is still good symmetry so as to suppress large corrections to the p parameter in the scale where 
the Lagrangian (3.33) is available. 
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The strongest unitarity condition is also given by Eq. (3.7) with the largest eigenvalue in the above. 

Cross sections of W^W£ — > hh and W^W^ — > W^W^ vanish because only a Higgs doublet has nonzero 
vacuum expectation value, and these interactions must include both of Z2 even and odd particles. 

In this model we have four coefficients Eqs. (3.38) and (3.39) to parametrize the dimension six differential 
operators. If we suppose that so = $l = vo = vl = 1, the eigenvalue c l 2 becomes the largest: c\ = 8. This value 
gives us the perturbative unitarity condition as Eq. (3.30). 

The unitarity bound, Eq. (3.7), can be interpreted as the perturbativity condition of couplings. For ex- 
ample, if so = sl = vo = and vl / 0, we get |C max | = 9vl/2. Then, the unitarity bound can be expressed 
as 

- 40T m vL 

VJ~ — • (3.46) 

3 §L 

In order to the preserve unitarity, the unitarity violation scale should be larger than the mass, m v i. As a result, 
we get the following condition: 




(3.47) 



If the model includes only one 3o vector, this requirement is stronger than the naive perturbativity condition, 
gL < 4tt. On the other hand, if it includes several 3o vectors, the unitarity bound also limits the number of 
them. 



4 Conclusion 

We have studied the perturbative unitarity for dimension six derivative interactions of the Higgs doublets. 
They are generated by kinetic terms in composite Higgs models or integrating out of heavy particles which 
interact with the Higgs doublets. Because of the latter case, derivative interactions appear even in models 
consisting of elementary Higgs doublets. 

We have firstly studied it in models including only one Higgs doublet. The strongest bounds for the 
unitarity are expressed by the largest eigenvalue of the matrix given by partial wave amplitudes of VBS 
processes. We focus on high energy region such that derivative interactions could dominate longitudinal 
mode contributions. Assuming that given derivative interactions respect the EWPM, only a combination of 
parameter, c H / f 2 , appear in the analysis. Therefore, the unitarity condition is express by the parameter as 
Eq. (2.6). We have applied it to the cross section of WlWl — > hh. 

We have calculated the bounds on explicit models: the minimal composite Higgs model; the littlest Higgs 
model. Structures of global symmetry are significantly different each other; SO(5)/SO(4) and SU (5)/50(5). 
However, given bounds are similar; c H = 1 and 1/2. Decay constants / are related to masses of the top like 
fermions in composite Higgs models. Therefore it is supposed that // Vc^ is larger than about 500 GeV, see 
e.g. Ref. [14], where the perturbative unitarity is violated above the region y/I > 2 TeV. Even in this case, it 
is difficult to obtain large cross sections enough to distinguish new physics contributions from the SM ones. 

Secondly, similar analyses have been performed in 2HDMs. A simple formula of the unitarity bound 
could not be obtained with parameters included in the effective Lagrangian (3.1) since the matrix of partial 
wave amplitudes are too complex to be diagonalized. Hence we have investigated the unitarity bound with 
explicit models: the bestest little Higgs model; the UV friendly little Higgs model; the inert doublet model. 
The first and the second ones are literally a kind of little Higgs model and the third one is a toy model including 
elementary Higgs doublets. 

In the first one, derivative interactions are generated by not only the kinetic term but also integrating out 
of a heavy scalar field. Constraints of the former one to the unitarity is similar to those given by lHDMs 
discussed in Sec. 2. Including the latter one, the largest eigenvalue depends on the scalar contribution. The 
unitarity bound can be stronger than the case not including it. 

In the second model, the unitarity condition is much severe comparing to the other models mentioned 
in this paper. This is because coefficients of derivative interactions are large in this model. These large 
coefficients produce large cross sections, so that the unitarity bounds of cross sections are also large. However, 
energy scales where the unitarity is broken down depend on the largest eigenvalue. Therefore the unitarity is 
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violated at low scale compared with that of other models. In this kind of model, large cross sections, enough 
to overcome the SM ones, can be achieved even with the large decay constant if people discard the naturalness 
criterion. On the other hand, assuming that masses of additional particles are near the decay constant, they are 
also near the scale of the unitarity violation. Therefore, study of vector resonance is probably required when 
people investigate, for instance, pair productions of these additional particles with vector boson collisions. 

In the last model which includes elementary Higgs doublet, VBSs of the SM particles are suppressed 
by the additional Z2 symmetry. In this kind of model, masses of heavy particles should be smaller than the 
unitarity violation scale. This condition means that couplings between Higgs doublets and heavy particles are 
much smaller than the strong coupling, 471, or the number of these particles are limited. 

As a conclusion, we have clarified the importance to study the unitarity bound when Higgs derivative 
interactions are investigated since the bound can be significantly lower than the naive cut off scale. 
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A Perturbative unitarity 

Amplitudes of elastic scattering satisfy the following relation for each partial wave: 

M l n =\{a,b)\M R n +iM l n ^ (A.l) 



X(a,b) = y/(l - (a + bf)(l - (a-b) 2 ), (A.2) 

where M^/M^ is the real/imaginary part of partial wave amplitudes, a and b are mass fractions between the 
mass of each particle and the center of mass energy, m a j, / ^Js, and partial waves are defined as below with the 
Legendre polynomials, P m {x): 



^(cosd) = 167r£(2« + l)M„/ 5 „(cos0), (A.3) 

.1 2 

dxP m (x)P n (x) = — — - 8 mn . (A.4) 
-1 2n+l 



1: 



Eq. (A.l) is the equation of circle whose radius is 1/2A and the center is (0, 1/2A). With the high energy 
limit where masses of produced particles can be neglected, the radius of the circle becomes the maximum. 
Therefore, the actual amplitudes are included by the maximal circle. Finally, partial wave amplitudes at least 
satisfy M* G [-1/2, 1/2] and G [0, 1]. If processes include the same particles in the final state, the bound 
becomes weaker as M„ <E [-1,1] andM^ £ [0,2]. 

Considering the unitarity bound for derivative interactions, with the massless limit, amplitudes of deriva- 
tive interactions can be written as 

- CJ+Qt (A.5) 



where s and t are the Mandelstam variables. The zeroth and the first modes of partial wave amplitudes appear: 



■Tf I IC-^r \P« + ^-P\ I, (A.6) 



m "=16^( c '-t)' <a ' 7) 

M ' =d?T (A - 8) 
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Eventually, the following conditions are obtained, respectively, 

i - l6n (A 9) 

If the condition, \C t \ < 3\C S \ < 2\Q\, is satisfied, the unitarity bound given by the first modes of the partial 
wave amplitude is stronger than that given by the zeroth modes. 

B Custodial symmetry of derivative interactions 

Dimension six derivative interactions of Higgs doublets naively violate the custodial symmetry. We study 
conditions which ensure the symmetry on derivative interactions for 2HDMs. In the following discussion, we 
refer results and notations in Ref. [15]. 

Derivative interactions are classified into three kinds of operators with combinations of indices: operators 
including unique indices are called type I, e.g. d (H\H\)d (hIH\); in type II, only one of four doublets has a 
different index, e.g. d(HjHi)d(HlH 2 ); the others belong to type III, e.g. d(HlH 2 )d(H%Hi). For tyP e 1 and 
II, current-current interactions, namely T operators, violate the custodial symmetry because they produce 
additional contributions to the mass of the Z boson. This is interpreted from a different viewpoint with the 
nonlinear representation. For type I, as studied in Ref. [16], operators belonging to T consist of the generator 
of the hypercharge i.e. the third generator of SU (2)r: 

(hT R3 dh)(hT R3 dh)= l -0 T nn , (B.l) 

where h is a real scalar multiplet which corresponds to the Higgs doublet. Since the generator violates 50(4) 
symmetry, the custodial symmetry cannot be preserved after the EWSB. In other word, operators which consist 
of 50(4) symmetric combinations of generators produce the custodial symmetric derivative interactions. It is 
also the case for operators of type II 12 . 

For type III derivative interactions, the situation is different, that is, certain combinations of SU (2)# violat- 
ing operators recover SU (2)r symmetry. In this type, the following operators produce real DOF of derivative 
interactions: 



r r , L(X '-pLcc 
1 11 1 22 


= 4(30i221 — Of 12 2 + Of 2 2l)> 


(B.2) 


tLol t'LCC 
1 12 1 12 


= 2(3^n22 + 30^ 12 + 30f 12 i + 0[ 122 - Of 221 ), 


(B.3) 


Ml i 22 


= 4 (3^1212 + 3O2121 + Of 122 - Of 212 - 2121 ), 


(B.4) 


1 12 1 12 


= 2(3^fi22 + 30^21 — Ofl22 + ^r212 + 0\ n \)i 


(B.5) 




= 2 (301122 — 30^> 12 - 30f 121 + Of 122 — 01221 )> 


(B.6) 


^12 ^12 


3 . £j jj^ H \ 
= 2(201122 — 01221 — 01212 ~ 0212l)) 


(B.7) 


1 n 1 22 


— 2^1122) 


(B.8) 


1 12 1 12 


= 2(01221 + 0r212 + 02121 )) 


(B.9) 


U 12 Un 


= 2 (01221 ~ 0f212 — 02121 )) 


(B.10) 



12 It is also true for imaginary DOF of derivative interactions. Any imaginary DOF of coefficients violate, in addition to the CP 
symmetry, the custodial symmetry. There are no relations so as to ensure the custodial symmetry for imaginary DOF, so that we 
discuss only real DOF here. 
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where Xf- := (hX^^dh) and h includes eight real scalar fields interpreted as two Higgs doublets. Generators 

Tfjjy S®fjj and Ifyj) are respectively (3, 1), (1,3), (3,3) and (1, 1) representations of SU (2)l x SU (2)r. 
Explicit forms of these matrices are given in Ref. [15]. If we naively follow the discussion for type I and 
II, coefficients of operators, (B.6), (B.8) and (B.9), apparently should vanish for SO(4) symmetry. However, 
we have found that operators preserving SU (2)r can be given by several combinations of operators violating 
SU(2)r. Since these operators are not linearly independent of each other, several relations are derived. In 
these relations, the above operators violating SU(2)r symmetry appear with a certain proportional relations: 

a 1212 : «1122 : «1212 = 1 : - 2 : !> ( B -H) 
where they are respectively coefficients of operators (B.6), (B.8) and (B.9). This condition is expressed as 

C[l22 + Cl221+Cl2l2=0, (B-12) 

3c[ 122 + cf 212 -cf 221 =0, (B.13) 

where c^y are defined as coefficients of O i - kl . The result is consistent with the custodial symmetric conditions 
shown in Ref. [15]. 

The above analysis is easily extended to models including N Higgs doublets. In the case, two other classes 
of derivative interactions should be defined: operators including three deferent indices are called type IV, e.g. 
d {Hj Hj)d (Hj Hk); the other operators whose indices are totally different each other are classified as type V, 
e.g. d(HjHj)d(H^.H[). With similar discussions given in the above, the following proportional relations are 
obtained: 



a ijik '■ a ujk '■ a ijik — 1 
a ijkl '■ a %kl '■ a Ijkl '■ a Ikjl '■ a Iljk = i 

S S Y Y Y 1 
a ikjl '■ a iljk : a ijkl '■ a ikjl '■ a iljk = 1 



-2:1, (B.14) 
— 1:1: —2 : 1, (B.15) 
1 : —1 : 1 : 1, (B.16) 



where 

®Jjkl ®fjkl respectively coefficients T^T^ and Sfj^Sff, and the first relation is for type IV and 

the others for type V. The following relations are induced for coefficients of derivative interactions: for type 
IV, 

4jk + cJ jik + cJj ki =0, (B.17) 

3 c Iijk + c ?jki - c fjik =0; (B . 1 8) 

for type V, 

c fjkl ~ c fjlk ~ c fkjl + c Mj + c Ujk ~ c flkj =0' (B-19) 

3(cjju + c fjik) ~ c fkji + c mj ~ c ujk + cflkj =0' (B.20) 

3(4;/ + c lij) ~ ^jki + cfjik + cfjk ~ cf kj =0, (B.21) 

3( c fij k + cfikj) + cfjkl - cfjik + cfji - cfij =0. (B.22) 

After imposing these conditions to ensure SO{4) symmetry on derivative interactions, remained DOF of them 
are shown in Tab. 2. 



C Unitarity matrices and other bounds 

Following matrices are given by zeroth modes of partial wave amplitudes for various VBS processes in 
2HDMs. Using the largest eigenvalue of them, the perturbative unitarity bound is obtained with Eq. (3.7). 



15 





with 


without 


I 


N 


2N 


II 


N(N-l) 


2N(N-\) 


III 


N(N- 1) 


3/V(/V- 1) 


IV 


N(N- 1)(JV — 2) 


3N(N- 1)(JV — 2) 


V 


/V(/V-l)(/V-2)(/V-3)/6 


/V(/V-l)(/V-2)(/V-3)/2 


Sum 


N 2 {N 2 + 5)/6 


N 2 (N 2 + 3)/2 



Table 2: Real DOF of dimension six derivative interactions on models including N Higgs doublets 
with/without SO(4) symmetry for each type. 

C.l Neutral two body states 

The matrix for partial wave amplitudes of neutral two body states are shown here. Initial and final states are 
given by eight states, namely, CjfCf , C^C 2 , C^Cf , C^C 2 ,NiN\, N r N\, N 2 n\ andA^vJ. If all of coefficients 
are, except for cf m , turned off, the matrix becomes the one in the case of 1HDM given in Eq. (2.5). 



/ 3c n 1 1 + c n 1 1 
2 

3c U12+ c U12 

2 

or , h 
££ni2+£hi2 

2 

3c 1221 ~ c f221 + 2c U22 
2 

M 
t llll 

c 1112 
Ji 

c 1112 

C fl22 
c l 11 1 

M 

c 1112 
„// 

c 1112 
„// 

c 1122 

-'Mill ' Mill 
2 

3c U12+ c U12 
2 

3c U 12+ C U 12 
2 

- 3C 122L M221 ' Z M12 

2 



££ni2_l£iii2 
2 



J M122 n 



Z M221 M122 



3c?,,,- 



-1112 



c 1221 



-1212 



M 
6 2221 



c 1 1 12 



-1221 



c 1212 



c 2221 

nT , H 
££lU22Mll2 
2 



3c U22+ 2c 122 



3c 1212+ c 1212 



'2221 



3c U12+ c U12 
2 



££1212211212 



J M122 ' Z M221 11 1122 
2 



JC 2221 ' c 2221 



-1112 



Ji 

c 1212 



-1221 



6 2221 



3c 1221 c f221+ 2c U22 



n T ,M 
- ,6 222l +<: 2221 



c 1 1 12 



-1212 



c 1221 



c 2221 

££lU2211112 
2 



3c 1717 _ t _c l91 



3c U22+ 2c f221 C U2 
2 



■ r ' c 2221 ' c 2221 



3c 2222+ c 2222 



-1122 



M 
c 2221 



-2221 



M 

c 2222 



122 



-2221 



c 2221 



C 2222 



V 7 " - 
-"- 1221 L 1221 



'2221 



'2221 



-2221 



-"-2222 ' c 2222 



(Cl) 



If we impose 50(4) symmetry on the above matrix and eliminate cf 122 and cf 212 with Eqs. (B.12) and (B.13), 
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the matrix is simplified as follows: 



2 

c 1 1 12 
2 

2 

-"-1221 c 1221 ' zc 1122 
2 

t llll 

c ni2 
„// 

c 1112 

C H22 
c 1 1 1 1 

c ni2 
„// 

c 1112 

c 1122 

*- 1 1 1 1 
2 

c 1 1 12 
2 

2 

JC 1221 L 122M ZC 112 

2 



C 1U2 
2 

c f221+ c 1212~ c U22 
2 



*- 1 1 12 
2 



JC I221 c 1221 



c 2221 
2 

L 1112 



6 1221 



Ji 
'1212 



6 2221 



-"-1221 c 1221 



c 1221+ c 1212 c fl22 
2 



c 2221 



^ C I221 c f221+ 2c U22 



c 2221 



c 1 1 12 



C 1221 



Ji 

c \2\2 



M 
-2221 



M 

*- 1 1 12 



^221 +^212 C U22 
2 



%r T -r" 
Jt-1221 c 1221 



Ji 
-1112 



c 1212 



C 1221 



c 222\ 



Ji 

c 1112 



C 1212 



Ji 
6 1221 



Ji 
6 2221 

c 1 1 12 
2 



-"-1221 c 1221 



c 1221+ c 4212 C U22 
2 



11 2221 



c 2222 



-1122 



Ji 
6 2221 



-2221 



Ji 

6 2222 



122 



C 2221 



Ji 
c 2221 



6 2222 

3c 1221 ~ c f221 + 2c U2 
2 



11 2222 



(C.2) 



C.2 Singly charged two body states 

The matrix for zeroth mode partial wave amplitudes of singly charged two body states is shown below. Initial 
and final states consist of the following four states: Cj l W 1 t ; C\n\; C^n\ and C%n\: 



f 3gfj { { +cf l { ] _^£m2+£uA2 3c U12+ c U12 3c 1221+ c 1221 \ 



3c U12+ c U12 3c U 22+^*1122 ^ c 1212+ c 1212 3c 2221+ i: 22; 



££iii2~n£iii2 



- >c 1221+ c 1221 



2 2 
Jc 1212 ' c 1212 - >C I122+ C 1122 



3c 2221+ c 2221 3C2221+ C 222I 



JC 2221 ' c 2221 



JC 2222 ' c 2222 



(C.3) 
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We explicitly impose SO(4) symmetry like the case of the neutral states: 
/ 



V- 



r" 
Mill 

2 

r" 
M112 

2 

11 1 1 12 
2 



3c [221+ c 1221 



_lui2 

2 

c 1221~ c f212~ c fl22 
2 



J M221 c 1221 



M 

M112 
2 



c 222l 



J< - 1221 c 1221 



<M221 c 1212 C U22 



c 2221 



3Cj22l+£l221 



c 2221 



c 2221 



,12222 



(C.4) 



C.3 Doubly charged two body states 

Entries of the following matrices are coefficients of VBS processes for doubly charged states: C^C^\ C^C^ 
and C^C^- If processes include the same particle states in their final states, i.e. the first and the third ones, 
the unitarity bound becomes weak as mentioned in App. A. The effect has been included in the following 
matrices: 



/- 


r T _ H 
Mill Mill 


C fll2 C fll2 


r T -r" 
M212 M212 


2 


2 




c 1112 M112 


3c 1221 + 3c U22+ c f221 +M122 


c 2221 t 2221 




2 


2 


2 


V- 


-T _ H 
c 1212 6 1212 


c 2221 c 2221 


c 2222 c 2222 


2 


2 



\ 



(C.5) 



We also impose SO (4) symmetry on the above one: 



Mill 
2 

r" 
M112 

2 



3c f22i M::i 



C1971 +4c^,, 



L 1112 

3c f221+ c {212+ c fl2 
2 



C 2221 



3c 



i ::l c f221+^ c 1212 



c 2221 



£2222 / 
9 / 



(C6) 
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